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ESTIMATING DRIFT PARAMETERS IN A FRACTIONAL ORNSTEIN 
UHLENBECK PROCESS WITH PERIODIC MEAN 

HEROLD DEHLING, BRICE FRANKE, AND JEANNETTE H.C. WOERNER 


Abstract. We construct a least squares estimator for the drift parameters of a fractional 
Ornstein Uhlenbeck process with periodic mean function and long range dependence. For 
this estimator we prove consistency and asymptotic normality. In contrast to the classical 
fractional Ornstein Uhlenbeck process without periodic mean function the rate of conver¬ 
gence is slower depending on the Hurst parameter iJ, namely . 
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1. Introduction 

The classical mean reverting Ornstein Uhlenbeck process is defined by the following sto¬ 
chastic differential equation 

dXt = —aXtdt + adBf, 

where a and a are positive real numbers and {Bt)t>o is a standard Brownian motion. Under 
Gaussian initial conditions the resulting process {Xt)t>o is a Gaussian process. This prop¬ 
erty and the strong mean reverting character of (At)t>o lead to a long history in various 
applications. Slight modihcations of the above SDE have been proposed to cover further 
features of observed data. For example we may add a deterministic drift term to the above 
stochastic differential equation; i.e.: 

dXt = {L{t) — aXt)dt + adBt- 

The resulting process is mean reverting around the time dependent center L{t). This may 
be used to model seasonalities or trends in the data. However, the process {Xt)t>Q is a 
semimartingale and thus the range of applications of the above models is limited to data 
which do not exhibit long range dependence. In situations with long range dependence but 
without mean reversion fractional Brownian motion (i?^)f>o is the classical model. Once 
the Hurst parameter H G (0,1) hxed it is characterized as the unique continuous zero mean 
Gaussian process with covariance 

V,[B«B«] = i (|t|“ + |s|“ - |i - sp") . 

Hence we see that for R = 1/2, we recover Brownian motion and for H > 1/2 we have long 
range dependence. One possibility of a representation of fractional Brownian motion is the 
moving average representation in Mandelbrot and Van Ness [5] of the form 

/ OO 

(ai((t - + a2((t - - {-sf_-"'‘'))dW,, 

■OO 

where = max(a:, 0), = max(—x, 0), Oi, 02 G M, R G (0,1), G is a normalizing constant 

and W denotes a two-sided Brownian motion. This is defined by Wt = if t > 0 and 
Wt = —if t < 0, where and W‘^ denote independent copies of W. 
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In view of combining the mean reverting character of the Ornstein Uhlenbeck model with 
the long memory property of the fractional Brownian motion we introduce for H > 1/2 a. 
generalized fractional Ornstein Uhlenbeck processes by the SDE 


dXt = {L{t) — aXtjdt + adB/^. 


( 1 ) 


First we have to specify in which sense we interpret integrals with respect to fractional 
Brownian motion. Unless in the case of Brownian motion, where due to the semimartingale 
setting the ltd integral is the canonical choice, this is not so clear for fractional Brownian 
motion with H > 1/2. We will now interpret our integrals as divergence integral (cf. e.g. [6]) 


since this type of integral allow us to derive a solution to (1) in an analogous form as in the 


Brownian case (cf. [2]). This would not be the case, if we interpret the integral as pathswise 
Riemann Stieltjes integral. Furthermore, in [4] it was shown that a least squares approach 
for the mean reverting parameter in a fractional Ornstein Uhlenbeck setting without periodic 
mean only leads to a consistent estimator when the integrals are interpreted as divergence 
integrals. 

In this paper we want to study the particular class of periodic drift functions L{t) = 
Y7i=i ) where the functions (pi{t)]i = 1, ...,p are bounded and periodic with the same 

period v and the real numbers = l,...,p are unknown parameters. Those parameters 
together with the positive real number a will be used to ht the model to the data. The 
parameter a > 0 will be assumed known, since for continuous observations it may be read 
of the sample paths immediately. Thus we have to estimate 6 := (/xi,...,/ip, a)* from an 
p + 1-dimensional parameter space 0 := x R’*'. Concerning the sampling scheme we 
assume that we observed the full path of the process (W)t>o ob a time interval T which is a 
multiple of the time period z/; i.e. T = nv for some n G N and let n —)■ oo. In the following 
the Lebesgue measure on the real line will be denoted by Without loss of generality we 
assume that the functions ipp^i = 1, ...,p are orthonormal in T^([0, z/], i.e. 



Moreover, we will assume that the functions pp i = 1, ...,p are bounded by a constant C > 0. 
The resulting SDE then takes the form 



with initial condition Xq = where is a random variable independent of the fractional 


Brownian motion 

The outline of the paper is the following. First we provide some properties of the model 
taking into account the special features of the periodic mean function. In section 3 we 
derive our estimator and in section 4 and 5 we prove consistency and asymptotic normality, 
respectively. 


2. Some preliminary facts on the model 

First of all we need some properties of our model which rely on the periodicity of the mean 
function. Now we derive an explicit solution for equation (1). 
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Proposition 2.1. The following stochastic process {Xt)t>o given by 

Xt = e-"* e^^L{s)ds - a ] t > 0 

is the unigue almost surely continuous solution of eguation (1) with initial condition Xq = .^o- 

Proof. This follows analogously to the Brownian case in [2] from the representation of the 
fractional Ornstein Uhlenbeck process with initial condition in [1], since the Ito formula 
for the divergence integral of fractional Brownian motion only differs in the term of the 
second derivatives, which do not contribute in the Ornstein-Uhlenbeck setting, cf. [6]. □ 

In the following we need a stationary solution of equation (1) in order to prove a necessary 
ergodic theorem. The stationary solution is given in an analogous way as moving average as 
in the classical Ornstein Uhlenbeck process. 

Proposition 2.2. The following stochastic process {Xt)t>o given by 

Xt := e-"* ( j e^^L{s)ds + a j c^^dBf 

\J —oo J —oo 

is an almost surely continuous solution of eguation (1). 

Proof. This follows from the representation of the stationary fractional Ornstein Uhlenbeck 
process in [1]. Note that we do not have any problem with the integral starting from —oo 
since the moving average representation dehnes the fractional Brownian motion on the whole 
real line. □ 


Next we show that for large t the difference between the two representations tends to zero. 
Proposition 2.3. As t ^ oo we obtain almost surely that \Xt — Xf —?■ 0. 

Proof. The explicit representations of Xt and Xt yield 

rO pO 


\Xt-Xt\ < e-“‘|eo|+e-“* 


e“U(s)ds 


+ ae 


—at 


e'^^dB 


H 

s 




where the right hand side tends to zero almost surely as t ^ oo. □ 

From this solution we may now construct a stationary and ergodic sequence of random 
variables which we need later for our limit theorems. 

Proposition 2.4. Assume that L is periodic with period 1, then the seguence ofC[0, l]-valued 
random variables 

Wk{s) := Xfe-i+s, 0<s<l,/ceN 

is stationary and ergodic. 

Proof. Since L is periodic, the function 

h{t) := e""* [ e^^L{s)ds 
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is also periodic on M. We have for any t E [0,1] that 


Wk{t) = e 




-k-l+t 


e“U(s)ds + 






= h(t) + 


>k-l 


= h{t) + ere 




s+k—l 


dB^ + (T y" f e^^dB 

z=-oo >^'-1 

V [ e^^dByi_^ 

_ Jo 


l=—oo 
0 


= h{t) + ae 


—at I „as J jaH 


+ a 


^ g-a(i+i-i) f e^^dBy^ 

Jo 


s+j+k— 2 ' 


]=-oo 


Thus, we have the almost sure representation 


Wfc(-) = M-) + Fo(W)+ ^ 


j=-oo 


with the functionals 


Fo : C[0,1] —)■ C[0,1]; a; i-A t i-A ere “ / e'^^doj{s) , 


F : ^[0,1] ^ C[0, l];a; ^ ere""* / e“*dn;(s) 

Jo 

and the (^[O, l]-valued random variable 

Yi := [s ^ 0 < s < 1] . 

The sequence of Gaussian random variables {Yi)i^z is stationary and ergodic. This implies 
that the sequence of G[0, l]-valued random variables {Wk)keN is stationary and ergodic. □ 


3. The estimator and its motivation 

In this section we want to motivate a particular estimator 6 for the parameter 6 through 
a least squares approach. The same kind of approach was used in [3] to derive an estimator 
for drift parameters in a Levy driven stochastic differential equations. First, we analyze the 
more general estimation problem of a p + 1-dimensional parameter vector 6 = {6i, 0p+i) 
in the stochastic differential equation 

dXt = ef{t,Xt)dt + adB^, 

where f(t, x) = (/i(t, x ),..., /p+i(t, x)y with suitable real valued functions fi(t, x);l < i < p. 
A discretization of the above equation on the time interval [0,T] yields for At := T/N and 
i = 

p+i 

(2) — Waz = y2 Xi^t)6jAt + a (i?^+i)Az “ • 
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The system of equations (2) has some resemblance to classical linear models with dependent 
noise. In accordance to the classical least squares approach for estimation in linear models 
we can try to minimize the least squares functional 

N / p+l \ 2 

(3) Q: (^1, ■■■, 9p+i) i-A I ~ XjAt ~ Xi^t)6jAt I . 

i=l \ j=l J 

It was discussed in [3] that the minimizer 6T,At of the above functional has the form 

^T,At = Qr,At-^T,Ai) 

with 

/ ... Eio/i(*At,EAi)/p(bAt,EAt)At \ 


Qtm = 
and 


V Eto/p(*At,X,Ai)/i(hAf,XiAi)At ... E!Io/p(*At,X,Ai)/p(hAt,X,At)At y 


N 


N 


PT,At ■ I ^ ^iAt){^[i+l)At At At) ; ^ /p+l (f At, A^Af ) i^{i+l)At AjAt) 


, i=l 


i=l 


This motivates the continuous time estimator 9x = Qj^Pt with 




Qt — 


CMt,Xt)fp+i{t,Xt)dt \ 


\ Jo" /p+i(t, A,)/i(t, A*)df ... /,+i(t,X0/p+i(t, A,)df 


and 


Pt : = 


h{t,Xt)dXt,..., U{t,x^)dXt 


In the special case of the fractional Ornstein Uhlenbeck process we have 6 = (/ii, ...,fip,a) 
and f{t,x) := ((pi, —x)^. This yields for T = nv the estimator 


(4) 

with 

and 


On := Q-^P. 


n n 


Pn-.= 


ipi{t)dXt ,/ ipp{t)dXt,- / XtdXt 


Qn ■ 


On 


G„-.= 


ir ••• ir 


= nulp, 


On : = 


ipi{t)Xtdt,..., / (pp{t)Xtdt 


where 
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and 


br, := / X?dt. 


Note that Ip denotes the p-dimensional unit matrix. 

First we deduce an explicit representation of the estimator 6n- For simplicity from now 
on we set z/ = 1. 

Proposition 3.1. We have 9n = 0 + aQ~^Rn with 


Rn ■ = 


(piit)dB^ / (pp{t)dB^ / XtdB 




Proof. Using equation (1), we obtain for the Fth component (1 < i < p) of the vector Pn 
the representation 


^Pi{t)dXt ='^ Hj / ipi{t)ipj{t)dt - a / ipi{t)Xtdt + a / ipi{t)dB 
Jo Jo Jo 


H 


and for the p + 1-th component 


XtdXt = V Pi / XtPj{t)dt -a X^dt + a XtdB^ 

Jo Jo Jo 


H 

t • 


This yields Pn = Qn9 + crRn from which together with 9 = Q^^^Pn proves the claim of the 
proposition. □ 


Furthermore, we can compute the matrix explicitly. 

Proposition 3.2. We obtain 

^ — 1 ^ f "Tn-^n \ 

~n\ -InK In )■ 

with 

/ I fn If"' 

An = (An,l, ..., An,p)* '■={- <Pl(t)dt, ..., - / ipp{t)Xtdt 

Jo ^ Jo 

and 

Proof. See [2] for the proof. □ 

Note that, since we will show in Proposition 4.2 that the limit of nQ~^ is well defined, 
especially lim„_,.oo7n = 7 > 0, this implies that for large enough n also nQ~^ is well dehned 
almost surely. 
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4. The consistency of the estimator 

In this section we prove that the estimator is consistent using the representation of 
Proposition 3.1. 

Theorem 1. For H G (1/2, 3/4) 9n defined in (4) converges in probability to 9 as n ^ oo. 

Proof. Since by Proposition 3.1 we have 9n = 9 + anQffi^Rn the statement of the theorem 
follows directly by the following two Propositions, Proposition 4.2 and Proposition 4.1. □ 

Proposition 4.1. For H G (1/2, 3/4) the sequence n~^Rn is bounded in Lfi. 

Proof. Assume that supj|<pj(t)| < C < oo, then by the isometry for fractional Brownian 
motion the variance of the first p components is 


Var { n ^ I ipi{t)dBfi ) = IE 


r*n 


.-H 


n “ / ipi{t)dB. 

0 

n pn 


pn pn 

= H{2F[ — l)n~^^ / / ipi{u)ipi{v)\u — v\^^~‘^dudv 

Jo Jo 


'0 JO 
- 2 H 


< 2C^H{2H / / {v - uY^-^dudv = Cfi 

Jo Jo 


In order to compute the variance of the last component, we write the solution of equation 
(1) in the form 

Xt = e + h{t) + Zt 


with 


m := e-“‘ Yi f*. f 

i=i -^0 


e°‘^ipi{s)ds 


and 


Z{t) := (Te-“‘ [ 

Jo 


—at I ^as j r)H 


e^^dBj. 


The variance of the last component then is 


Var I n-^ I XtdBf ] = 


f*n \ 2 

XtdB^ 


'0 


= n“2^1E 


(e + h(f) + dB 


— El + E2 + E/3 + E/4 + E/5 + E/g, 
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with 


El := 

E 2 := 2 n- 2 ^E 


e-^%dB! 


'0 
pn 

Jo 
1 /•! 


= E 


1 rl 



e-»J^+J\u-v\^^-^dudv, 


0 ^0 


e-'^^iodB^ / h{t)dB^ 


H 


= 2 E[^o] [ [ ^ “"“/i(n-u)|M — “^dudv, 

Jo Jo 


E 3 := 2n-2^E 
E4 := 

Eg := 2 n- 2 ^E 


Jo 


e-'^^^odB^J / ZtdB, 


■H 


h{t)dB.^ 


H 


'0 

rn 


1 pi 



h{nu)h{nv)\u — ‘^dudv, 


0 Jo 


I h{r)dB^ / Z{t)dB^ 
Jo Jo 


H 


Ee := n“^^E 

( r z,dB«) 

= n-^^cj^E 



\Jo J 


\Jo Jo J 


and 


We now have to discuss the boundedness of each of those terms separately. The term Ei 
is obviously bounded as n —)■ cx). Since the function t i-A h(t) is bounded, it follows also 
that the terms E 2 and E 4 stay bounded. E 3 and E^ are zero by the isometry property for 
multiple Wiener integrals of different order (cf. e.g. [7]). Eq also stays bounded as n ^ 00 , 

since we know by [4] that ZtdB^Y is convergent for H G (1/2, 3/4). □ 

Proposition 4.2. As n ^ 00 we obtain that nQ~^ converges almost surely to 

Ip + jAA^ - 7 A 


C : = 


— 7 A* 


7 


where 


and 


r*l pi N ^ 

A = (Ai,..., Ap)* ^ (pi{t)h{t)dt, J (pp{t)h{t)dt 


7 : = 


[ h\t)dt + a^a-‘^^Hr{2H) 

JO i=i 


with h{t) := e /!oo 

Proof. We use the following notation 

w = hit) + Zi 


h(t) := e / e“V*('S)ds 

i=i J-°° 


with 
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being periodic with period 1 and 


Zt := ae 


—at 




We investigate the limit behavionr of the different entries in n ^Qn separately. 
The ergodic theorem, Proposition 2.4 and Proposition 2.3 yield 


lim kn,i = 

n^oo 


= IE 


rn 2 r'^ ^ ^ 

lim — / Xt(pi(t)dt = lim — / Xt^i(t)dt = lim — \ / Xt^i(t)dt 

n^n n Jq n^n n Jq n^n Tl ^ J 

t 

dt, 


Xt^pi{t)dt 


L-'o 


= / h{t)^pi{t)dt + / IE 

Jo Jo 


uo 


where the second term in the last expression is zero by the properties of Wiener integrals. 


Note that the seqnence of C[0,1] valned random variables 


s I —y Z, 


k-\-s 


and ergodic, hence by the ergodic theorem we have as n —)■ cxd 


; G Z is stationary 


- [ Ztdt = -'^ [ 


IE 


n 


Zfdt 


uo 


Moreover, we have 

1 r 


Zf — ZA dt 


n Jo 

Those two facts now imply that 


<T i 

<- / e 
n Jo 


—at 


e^^dB^ 


dt 


— IE[Zo] — 0. 


0 as n —)■ cxo. 


lim 

n—^oo 


Zfdt 


= 0 . 


It then follows from the bonndedness of the fnnction t ha h{t) that 

1 


lim snp 

n^oo 


1 r 

— I Xtdt 


n 


’0 


= lim snp 

n—>-oo 


n 


(e + h{t) + Zj) dt 


< oo 


From Proposition 2.3 we also have 


lim snp 

n—>-cx) 


Xfdt 


n 


< oo. 


Those two inequalities together with Proposition 2.3 yield 


- [ X^dt-- [ X^dt 

_ 

- f (Xt + Xt) (xt - Xt) dt 

n Jo n Jo 


n Jo \ /V J 


0 as n —)■ oo. 


Now nsing the ergodic theorem for h and [4] for Z we obtain 


1 1 
lim - / XUt = lim - ' 

n—>-00 n _/g n->-oo Tl 


Xfdt 


>0 




= / JA{t)dt + a^a-^^Hr{2H), 

Jo 

noting that the mixed term is zero dne to the properties of Wiener integrals. 
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Thus the general expression for 7 „ in Proposition 3.2 yields 

-1 

/ 1 ^ ^ ^ 

lim = 

n—¥oo n—>-OD 


hm i rXfdt-^Kl 

pt V 

Jo .-I 


-1 


Since the functions 997 i = 1, ■■■,p are orthonormal in T^[0,1] we can use the Bessel inequality 


p ^ ~ \ ^ z”! 

= ipi{t)h{t)dt] < 

i=i i=i / >^0 


h‘^{t)dt 


to see that the above limit is well dehned and hnite, which completes our proof. 


□ 


5. The asymptotic normality of the estimator 


In this section we prove asymptotic normality of our estimator which may be reduced to 
a limit theorem for dependent normally distributed random variables. 

Theorem 2. For H G (1/2, 3/4) we obtain for idn defined by (f) 

- i9) A AA(0, a^CSoC) 


with 

where 
G := 


Sn := 


G —a 
* b 


—a 


a* := I an 


10 JO r-i' 

fo' /o 

r/' 


IWI^- 

'^dsdt ... 

Fi{s)pp{t)\t- 

- ‘^dsdt ^ 


s\^^-^dsdt ... 

Id Id Fp(^)Fp(i)I- 

-s\^^-^dsdt ) 

3)h{t)\t- 

- ‘^dsdt, f f Fp{s)h{t)\t— 

Jo Jo 

s\‘^^~^dsdt^ , 

b := an 

[ [ h{s)h{t)\t 
Jo Jo 

— ‘^dsdt, 



an = H{2H- 

1), 


h{t) 

P pt 

i=l 

e^^Pi{s)ds 



and C is defined in Proposition (4-2). 


Let us discuss the differences to the Brownian case (cf. [2]) and the fractional Ornstein 
Uhlenbeck case without periodic mean function (cf. [4]) before proceeding with the proof. 
The rate of convergence is slower than in the Brownian case. Furthermore, it is 

also slower than the rate for the mean reverting parameter in a fractional Ornstein 

Uhlenbeck setting with L = 0. This is due to the special structure of our drift coefficient, 
which in our setting also dominates the component of a leading to a slower rate even for 
a and a different entry in the covariance matrix. Furthermore, unless in the Brownian 
case So C~^. This is due to the isometry formula for fractional Brownian motion with 
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H > 1/2, which is not simply derived from the scalar product in but from the scalar 
product in a larger Hilbert space H. Namely for a fixed time interval [0,T] T-L is dehned as 
the closure of the set of real valued step functions on [0, T] with respect to the scalar product 

< l[o,s] >'H= B^). 

Proof. By the representation 

= aQ-^Rn 

and the almost sure convergence of nQ~^ —)■ C it is sufficient to prove that as n —)> oo 

t 


n 


-H 






iPp{t)dBl^, —n 


.-H 


XtdB, 


H 


V 


AA(0,Sc 


As in the proof of Proposition 4.2 we may replace Xt by Xt with the representation Xt = 
Zt + h{t) and may deduce that Zt does not contribute to the covariance matrix. Namely 
the contributions to the off-diagonal elements in d and the mixed term of b are zero by the 
isometry formula for multiple Wiener integrals of different order (cf. e.g. [7]). Furthermore, 
Var(n“^ ZfdB^) —)■ 0 as n —)■ cxo, since we know by [4] that ^Var( J” ZtdB[^) is convergent 
and 2H > 1 for 1/2 < 77 < 3/4. 

Hence it is sufficient to show that for the 1-periodic functions (1 < i < p) and h as 
n —)■ oo 

t 


n 


-H 


Pi{t)dBf,...,n 


.-H 


^Pp{t)dBf, —n 


.-H 


h{t)dB. 


H 


V 


AA(0,Eo). 


□ 


This is an immediate consequence of the following Proposition 5.1. 

Proposition 5.1. Let /fc (1 < fc < m) be periodic real valued functions with period 1, then 
for H > 1/2 and n ^ oo 

t 


n 


-H 


fi{t)dBf,...,n 


.-H 


fm{t)dB, 


V 


Af f 0, 77(277 - 1) mfjis)\t - s\^^-^dsdt^ 


i 


Proof. Since fk is periodic with period 1, we may write for 1 < 7 < m 


n 


.-H 


pn pi ^ 

Jo Ji -1 


with 


and 


1 pi 


Y/^ ^Afi 0,77(277- 1) 



fk{t)fk{s)\t - s 


0 JO 


^^-^dsdtj 


Cov(y,^y') = pH{\t-j\)H{2H-i) 


1 



0 Jo 

1 


n 


2H-2 tt2 


H\2H-iy 



fk{t)fi{s)\t- ^dsdt 
fk{t)fi{s)\t - sy^~‘^dsdt 


0 Jo 


for 1 < f, j < n and 1 < 7, / < m, since 

1 
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Hence the sequences satisfy the conditions of Theorem 7.2.11 in [8], which implies as 

n —)■ cxo 

n-^ ^ r/ A AA (^0, H{2H - 1) fk{t)Ms)\t - sf^-^dsd?j . 

Finally the Cramer-Wold device together with a similar argument for the covariance terms 
lead to the desired result. □ 
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